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We investigate the phase transitions in two- legs ladder systems in the incommensurate phase, for 
which the gap is destroyed by a magnetic field (h c i < h) and the ladder is not yet totally saturated 
(h < h C 2). We compute quantitatively the correlation functions as a function of the magnetic field 
for an isolated strong coupling ladder J± 3> Jy and use it to study the phase transition occuring in 
a three dimensional array of antiferromagnetically coupled ladders. The three dimensional ordering 
is in the universality class of Bose condensation of hard core bosons. We compute the critical 
temperature T c (h) as well as various physical quantities such as the NMR relaxations rate. T c has an 
unusual camel-like shape with a local minimum at h = (h c i +h C 2)/2 and behaves as T c ~ {h — h c ±) 2 ^ 
for h ~ h c i. We discuss the experimental consequences for compounds such as Cu2(C5Hi2N2)2CLi 



I. INTRODUCTION 



There has been recently considerable interests on spin 
ladder materials. These systems, quite remarkably, have 
a gap in the spin excitation spectrum for an even num- 
ber of legs and no gap for an odd number. This_phc- 
nomenon, reminiscent of the Haldane conjectureHa has 
been explored p. gxeat details both theoreticallyQliJ and 
exp er iment allylij E3 . 

In ladders, contrarily to the case of spin S chains, the 
gap and the dispersion in the ladder are controlled by 
two different energy scales, namely the transverse Jj_ 
and longitudinal J|| exchanges. The ladders are thus 
prime candidates to study quantum phase transitions 
where the spin gap is destroyed by application of a mag- 
netic field. Because of this separation of energy scales 
between the gap and the exchange, even when the gap 
is destroyed quantum effects are still crucial. The lad- 
ders thus offers the possibility of an extremely rich quan- 
tum behavior, unsuspected in more conventional spin sys- 
tems. Such quantum, pliase transitions were indeed stud- 
ied experimentallyoO'ES. On the theoretical side^-they 
were investigated using a bosonization techniqueO for 
a single ladder. Close to the critical point where the 
gap vanished, the spin-spin correlation functions were 
found to diverge with a universal exponent, leading to 
a divergent NMR relaxation rate 1/Xi ~ T" 1 / 2 , in good 
agreement with the experimental findings. Between the 
critical field h c \ where the gap was destroyed and the 
saturation field h C 2, the ladder had incommensurate spin 
spin correlation function with a quite distinctive spec- 
trum compared to single chain systems. These results 
were confirmed and extended—in subsequent analytical 
and numerical calculationsc3~E3. 

Due to the gaped nature of the excitations for a sin- 
gle ladder when h < h c \ a weak interladder coupling is 
inefficient and the single ladder approximation is nearly 



exact. This is clearly different in the incommensurate 
phase h c \ < h < h C 2, and the question of the coupling 
of ladders becomes much more crucial. Quite generally 
interladder coupling can lead now to a three dimensional 
ordered phase. This is the case for example for the com- 
pound Cu2(C5Hi2N2)2Cl4 which has an experimentally 
accessible gap of A ~ 11K. Specific heat measurements 
have revealed the existence of a transition at fin it e t |epr 
pcrature, the nature of which is still controversialli§EJ~E!3. 
It is thus a challenge, both from a theoretical point of 
view and in view of application to experiments, to un- 
derstand how three dimensional ordering can occur in 
ladder systems. 

We investigate the nature and physical properties of 
such transition by looking at antiferromagnetically cou- 
pled ladders. Because of the peculiar nature of the excita- 
tion spectrum in ladders, this transition is different from 
the one occuring in more conventional spin materials. 

The plan of the paper is as follows. In section [il] wc 
define the model for coupled ladders. For simplicity we 
confine ourselves to the case of strongly coupled ladders 
J± 3> Jy . In section III we examine the single ladj 
this limit, using a mapping on a single spin chaired 
We compute quantitatively the correlation functions as 
a function of the magnetic field. The weak and strong 
coupling limits give an identical structure for the corre- 
lations functions and we recover the universal exponents 
and spectrum for the spin-spin correlation functions de- 
rived in Ref. [H| The three dimensionally coupled lad- 
ders are described in section IV. The three dimensional 
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ordering is in the universality class of Bose condensa- 
tion. We compute the critical temperature T c as well as 
various physical quantities such as the NMR relaxations 
rate. T c has a camel-like shape with a local minimum 
at h = (h c i + h C 2) and behaves as T c ~ (h — h c i) 2 ^ for 
h ~ h c i. We discuss the experimental consequences for 
compounds such as Cu2(C5Hi2N 2 )2Cl4. Conclusions can 
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be found in section [v] and some technical details are left 
for the Appendix. 



II. THE MODEL 

We consider the two legs ladders shown in Figure [j]. 



in 

FIG. 1. The two leg ladder system considered in this pa- 
per. An interladder coupling J' couples the ladder in a three 
dimensional way. 

For the moment we consider a single ladder and thus 
take J' = 0. The ladder Hamiltonian is given by 

H = J|| ^2 Si,t ' + J-L ^,1 ' ~ ^ 

(1) 

where I = 1, 2 denote the two legs of the ladder, and /i 
the applied magnetic field. 

The case when the rung coupling J± is much smaller 
than the interaction J\\ along the ladder has been studied 
by a variety of techniques both in thfirabsence o&T or 
in the presence of a magnetic fieldl!2E3 We concentrate 
here on the opposite limit J± ^> Ju . In that case the 
der can be mapped onto a single spin 1/2 chainE3 ! t 
and we recall the mapping here for completeness. Indeed 
an individual rung may be in a singlet or a triplet state. 
Applying a magnetic field brings one component of the 
triplet closer to the singlet ground state such that for a 
strong enough magnetic field we have a situation when 
singlet and m = — 1 component of triplet create a new 
effective spin-1/2. It is thus possible if J± ^> Jm to retain 
only these two states for all the magnetic field range be- 
tween h c \ when the gap is broken to h C 2 when the ladder 
is completely magnetized. 

One can easily project the original Heisenberg Hamil- 
tonian (m) on the new singlet-triplet subspace 



ll> = ^[| U)-|IT>] 
If) = I TT) 



(2) 



This leads to the definition of the effective spin 1/2 op- 
erators 



~~ V2 
Sl. 2 = \[I + 2~S z ] 



(3) 
(4) 



When expressed in term of the effective spin operators 
(If); the original Hamiltonian (]l|) becomes 

H c b = J|| y~^[5f Sf +1 + Sf Sf +1 + -xSfSf^i] 

i 

-hJ^Sf+C (5) 



where C = {— 4 L + lf — ^)Lisa simple energy shift and 
the system is in an effective magnetic field 



Jm 

h = h — J \ 

2 



(6) 



The Hamiltonian (|5|) describes a single spin 1/2 chain 
with a fixed XY anisotropy of 1/2 in an effective mag- 
netic field. In the following we denote with a tilde, the 
magnetic field h, and the magnetization rh of the effec- 
tive spin 1/2 chain. The gaped phase h < h c \ for the lad- 
der corresponds to the negatively saturated magnetized 
phase for the effective spin chain, whereas the massless 
phase for the ladder corresponds to the finite magnetiza- 
tion phase for the effective spin 1/2 chainEM The field h C 2 
where the ladder is totally magnetized correspond to the 
fully magnetized phase for the effective spin 1/2 chain. 
It is easy to check that 
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(7) 



III. SINGLE LADDER 

Before taking into account interladder interactions let 
us first recall some important consequences of such a 
mapping for the single ladder. In the process we give 
a more quantitative calculation for the correlation func- 
tions as a function of the magnetic field. The results of 
this section will be used to study the interladder cou- 
pling. We focus here on the massless phase h c \ < h < 
h C 2- To conveniently derive the low energy properties 
of the effective spin 1/2 chain we use the by now well 
known bosonization technique. We refer the reader to 
Rcf. ^7] |29|,[l9] for details and just recall here the main 
steps. r-n-. 

We first use the Jordan- Wigner transformationE§~Ell 
which essentially maps the spin problem onto a prob- 
lem of interacting fermions on a lattice. For the spin 
1/2 system considered here, the corresponding fermionic 
problem has Fermi momentum kp — -| if h = 0. Finite 
magnetic field corresponds to a chemical potential for the 
fermions. We then perform a linearization around the 
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free Fermi points given by ±kp , to obtain an effective 
low energy continuum fermionic theory and then express 
the fermion operators in term of bosonic ones related to 
the fermion density fluctuations using the standard dic- 
tionary of Abelian bosonization. 



S+(x) = 



9(a) 



S z {x) = — d x (, 

7T 



l ^ + cos2<p(x)] 
lS * cos2</>(x) 



(8) 



Where S + (x) — ^s=, S z (x) — ^ for x = na, a being 
the distance between two nearest neighbors sites along 
the chain. From now on we take the lattice spacing 
a = 1 and measure all distance in units of a. The 
phase 4> is related to the average density of fermions 
(or equivalently to the uniform spin density along z) by 
S z (x) — — ~d x <f), whereas 9 is connected to the conjugate 
momentum II of 4> (such that [4>(x),H(x')] = i8(x — x')) 
by 9{x) — f* dyll(y). In a very crude sense <j>, 6 can 
be viewed as the polar angles of a spin. The low energy 
properties of the Hamiltonian (JsJ) can be totally described 
in terms of the boson Hamiltonian 



H 



dx 
2^ 



U K(7Tll) 2 + -{d x <t>) 2 
K 



(9) 



where 
tion 



has been shifted to absorb the finite magnetiza- 



<f> — > cf> + 2mx 



(10) 



The only two parameters controlling the low energy prop- 
erties are the "spin wave" velocity u and a number K 
called the Luttinger liquid exponent. Both are known 
exactly for the spin 1/2 chained. For h — analytic ex- 
pressions are known 



Jz/Jx,y = - COS7T/3 : 

l/K = 2/3 2 
1 



(11) 



j _ -simMl-/? 2 ))^ 



Thus K = 1/2 for an isotropic Heisenberg chain with 
h = whereas K = 1 for the pure XY one. For the 
Hamiltonian this leads to 



K = 3/4 



3V3 J\\ 



(12) 



At finite magnetic field K and u can be obtained by inte- 
gration of the Bethe ansatz equations and are shown on 
Figure|| for the specific case of the XY anisotropy 1/2. 
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FIG. 2. Magnetic field h and Luttinger liquid parameters 
u and K for an XY anisotropy of 1/2 plotted as a function 
of the magnetization rh (only positive values are shown, the 
parameters being symmetric with rh — > —rh). rh — 1/2 is the 
saturated chain. K = 3/4 for zero magnetic field, whereas 
K — > 1 and u — > close to saturation since the excitations 
above the ground state become very diluted, u and h are in 
units of Jy/2. 

Close to h c i or h C 2 the number of excitations compared 
to the fully polarized ground state becomes very small (in 
the fermionic language one is close to an empty or a full 
band) , and thus K take the value for noninteracting par- 
ticles K — > 1 regardless of the strength of the original 
interaction J z /J xy . 

Since the free boson theory given by ([)]) is trivially 
solvable, it is straightforward to calculate the asymptotic 
decay of the dynamic correlation functions, which are just 
the ones of a spin 1 /2 chain. Using (||) , one gets for T = 
(for more details see e.g. Rcf. |l^) 



i / r 2K 



{S z (x, t)S z (0, 0)) = m 2 + Ci^ + C* 2 cos(tt(1 - 2m)x) 
(S + (x,t)S~(0,0)) = C 3 cos(2nmx) l-\ + C 4 cos(ttx) ( - J (13) 



where r = x 2 + (ut) 2 and Ci are constants on which we focus later in this section. When expressed in term of the 
true magnetization 2m = 1 + 2m and the original spin operators of the ladder using (§) this gives (e.g. for rung 1) 



1 fV 2K 
(S z (x,t)S z 1 (0,0)) = — + — +cos(27rmx) ( - 

2K+l/{2K) A\l/(2ff) 

(^(x,^r(0,0))=cos( 7 r(l-2m).x)(^j +cos(ttx)(-\ (14) 



(where we have dropped the constants C for simplic- 
ity). Equ. (|l4|) presents some remarkable features. First, 
as already pointed out in Ref. |l^, low energy modes ap- 
pear only close to q ~ in the S z correlation function 
or close to q ~ n for the transverse one. The q ~ tt 
(for S z ) or q ~ (for $ ) mode remain massive. This 
is in marked contrast to what would happen for a gaped 
(e.g. dimerized or frustrated) single chain system (in 
weak coupling) where both the q ~ and g ~ 7r would 
become massless when /i > h c \ . A summary of the mass- 
less and massive modes is shown on Figure || 
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(a) (b) 
FIG. 3. Schematic picture of the field dependent disper- 
sion as seen by the S Z S Z (a) and S + S~ (b) correlations for 
fields close to h c i. Only the dominant singularities are shown. 
In marked contrast with a single chain there is no massless 
excitations close to q ~ tt (resp. q ~ 0) due to the presence of 
an antisymmetric massive mode, unaffected by the magnetic 
fieldN. 



Such prediction for the ladder correlations should be 
testable in neutron or resonance experiments. Close to 
h c i (or h C 2) K — » 1 and one recovers the universal expo- 
nent for the decay of the correlation functions predicted 
in Ref. [l9| The weak coupling approximation only al- 
lowed for a qualitative calculation of the exponents far 
from h c i and h C 2- For the strong coupling case one can 
get the full magnetic field dependence as shown in Fig- 
ure ^. As shown in Appendix the whole asymptotic 
structure of the correlation function is independent of the 
strength of the coupling J± vs Jy provided of course that 
the correct Luttinger liquid exponent are used. 

Let us now compute quantitatively the correlation 



functions. We only focus here on the massless modes. 
The prefactors in the correlation functions of the spin 
operators in the XXZ model have been computed by 
Lukyanov and ZamolodchikovHEj. Thus for example for 
the transverse staggered magnetization we have 

{-l) n S + (x = no, r) = [Fp/^ia/uf^e* 6 ^ (15) 
where the expression for the prefactor readsB 

/3 2 



F 



1 



2(1 - /3 2 ) 2 



T{-B. 



2-2/3 2 - 



exp 



2V^r( 



dt 



2-2/3^ 



(16) 



sinh(/3 2 t) 



t Vsinh?jcosh[(l - /3 2 )t] 



In the vicinity of the value of interest here /3 2 
analytic calculation is possible: 



2/3 an 



Fa = 



2tt 2 / 3 



£(2/3) 

r(i/3) 



= -(/3 2 -2/3)[7+ln(37r/16)+ 7 r/v / 3] 



i^(2/3)exp[-2.173(/? 2 - 2/3)] 



(17) 



For more general values of K and m the value of F is 
shown on Figure |] 
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FIG. 4. Prefactor F as a function of the Luttinger pa- 
rameter K. The divergence close to K = 1/2 is due to the 
appearance of logarithmiC|-eorrections in the correlations func- 
tions for the isotropic caseEl For the specific case of the chain 
with XY anisotropy 1/2, the prefactor is shown as a function 
of the magnetization m. Note the relatively weak dependence 
on the whole interval due to the finite XY anisotropy of the 
effective spin chain. 

Using (|l^) one gets for the staggered part of the trans- 
verse spin-spin correlation function of the physical ladder 
spins as a function of space, time and temperature 



(^(no,r)S-(0,0)> = H(£) 



nT 



sin7rT(r + ix/v) 



(- 
(18) 



The term (—1)^ indicates that due to the relation be- 
tween physical and effective spins (||), the physical cor- 
relation function is singular at the wave vector (7r,7r). 
Let us insist that-fche position of the singularity does not 
change with fieldl!3 (see (14)). 

Assuming that the local susceptibility is dominated by 
the contribution from the transverse part of the stag- 
gered susceptibility, we get the following expression for 
the NMR relaxation rate: 



(19) 



(up to the hyperfine coupling constants). This leads, 
when the Fourier transform and analytical continuation 
of (fl8|) is performed (see also (|27|)), to the relaxation rate 



1 



F, 



oc ^r 2 (/?72)r(i- 1 ) 



T 



At (3 2 = 2/3 (i.e. for h 
gives for T\ 



(2ttT) 2 - 
or h = (h c i - 



(20) 



h c2 )/2) this 



(21) 



jn/3 



/3 3\/3r(2/3) 
16V2tt 



0.1 



T l/3 



1/3 



(22) 



Close to h c i or h c2 , K ^ 1 (thus (3 2 1/2) and © 
gives back the universaO divergence of the relaxation 
time 



T * jn/2 



(23) 



Away from the critical field the exponent increases 
weakly to -1/3 at h = 0. (see also Ref |l|,||). The full 
magnetic field dependence can be obtained from Figure |^ 
(using P 2 = 1/2K and ©). 

Although, as shown in Appendix the correlation 
functions for the strong and weak coupling ladder are 
smoothly connected, a very interesting question is how 
the Luttinger liquid parameter varies with the field when 
going from weak to strong coupling. This is not trivial 
since for the weak coupling ladder when A -^i h -^i 
one recovers essentially the Luttinger liquid exponent of 
a single chainEj. For an isotropic system this is K = 1/2, 
i.e. the same value than the universal one close to h c \. 
If there is XY anisotropy the parameter K increases and 
thus K > 1/2. On the other hand for the strongly cou- 
pled ladder regardless of the XY anisotropy (provided 
J z > 0) the luttinger liquid parameter decreases with the 
field (see (|,0) and Figure §), giving the very different 
field dependence shown in Figure |^. 
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FIG. 5. Qualitative variation of the Luttinger liquid pa- 
rameter If as a function of the magnetic field both for the 
weak and strong coupling ladder. Full line is the weak cou- 
pling ladder J± <C Jy , and dashed line for the strong cou- 
pling one J± ^> Jii . Different curves correspond to different 
XY anisotropy from the isotropic case to the XY limit. A 
is the gap for the weak coupling case. See also Figure Q for 
the exact result for the strong coupling ladder (due to the 
difference of definition (see text) between A" wea k and K s t r0 n S 

We plot A" st rong /2). 

It would be extremely interesting to have more quan- 
titative estimates for the behavior shown in Figure |^. 
An amusing possibility would be to have, for a certain 
strength of coupling and a given anisotropy a luttinger 
liquid parameter totally field independent. 
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IV. COUPLED LADDERS 

In order to describe realistic compounds we now take 
into account an interladder interaction of the form shown 
in Figure |l| and given by the Hamiltonian 

H 3D = ^ ^ladder + J ' X X S ' l - a ^ ' Si ^' 2 
at (a, (3) i 



where (a, B) denotes pairs of nearest neighbors ladders. 
It is easy to see from Figure [j] that a spin on leg 1 of 
one ladder can only interact with the spin on leg 2 of the 
neighboring ladder and vice versa. 

Since the interladder coupling is very weak it is again 
legitimate to map the problem to an effective spin 1/2 
problem. The coupled ladder system thus reduces to a 
problem of spin 1/2 chains coupled by the interaction 



^coupling = "J E + h ' C -] + l£^ + T^ (25) 

(a,p) (a, /3) a 



where z is the coordination number. Because a spin 
on leg 1 can only be coupled to a spin on leg 2 by J' 
this leads to a ferromagnetic coupling for the XY part 
of the interchain coupling although the original interlad- 
der coupling is antiferromagnetic. There is also a trivial 
redefinition of the effective magnetic field by the interlad- 
der coupling. Although the problem of coupled ladders 
now look identical to the one of three dimensionally cou- 
pled spin 1/2 chain, the. physics will be quite different 
from the standard caseE3 of isotropic coupled spin 1/2 
chains. Indeed, as we will see below, the XY anisotropy 
of the effective spin 1/2 chain inherent in the ladder prob- 
lem, plays a crucial role. The treatment of (^5|) depends 
crucially on what is the characteristic energy scale of the 
interladder coupling when compared to what happens for 
a single ladder. 



A. High density limit 

If one is far enough from h c x and /i C 2, the interladder 
coupling will be small compared to the characteristic en- 
ergies (Fermi energy for the associated spinless fermion 
problem) of the single chain. It is then possible to treat 
the interladder coupling in the mean field approximation 
while keeping the full single ladder physics. Since the 
single chain correlation functions along z decays faster 
than the one in the XY plane (see e.g. (|l4]), three di- 
mensional order will occur first in the XY plane. It is 
thus possible to neglect in fl25| ) the interchain S Z S Z cou- 
pling, and to retain only the XY part. Note that in that 
case it is not important whether the interchain coupling 
is ferro or antiferro since one can go from one to the 
other by making the gauge transformation S x — > — S x , 
S v — » —S y , S z — > S z on alternating chains. The coupling 
is just a spin flip term for the spin which in a bosonic 
representation for the chain is just a hoping term for the 
bosons (see below). Another way of viewing it is as a 



Josephson coupling between the phases 9 of the spins 
on different chains. This shows that the transition is 
a normal-superfluid type transition or alternatively is a 
Bose-condensation transition for the hard core bosons as- 
sociated with the effective spins. In a pictorial level this 
says that the orientation of the spins of the ladder in the 
XY plane want to lock in the same direction as shown in 
Figure ^. Since we are in d = 3 (i.e. .above the critical 
dimension for the quantum transitior£3) the exponents 
are the mean-field ones (y = 1/2 and £ = 1). 
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FIG. 6. Cartoon of the three dimensional transition in 
ladder systems. The direction of the spins in the XY plane 
tends to lock together between different chains leading to a 
planar antiferromagnet. Note that the triplet states are in 
fact delocalized on the ladder and that in the XY plane the 
spins remains modulated at q = 7r in the ladder direction. We 
have represented the singlets as shaded boxes. This transition 
is in the universality class of normal-superfluid transition or 
Bose-condensation of hard core bosons. 

To compute the transition one can use a standard 
mean-field approach. The transition temperature is given 
by 



1 

J 7 



X±(<Z = 0,u; = 0,r) 



(26) 



where x± is the single chain transverse staggered sus- 
ceptibility (wavevector q is to be counted relative to 
7r). For thft-Tuttinger liquid the susceptibility can be 
computedE§El 



X±(u,q) = — g-- 



B 2 B 2 
^—)B(^— 
2 4 



is-, 1 



2 ; 



l-/3 2 /2j 



(27) 



G 



where 



vq 



4ttT 



(28) 



and B(x,y) = T(x)T<j 
Solving © with (| 



/r(x + y). 

) gives the critical temperature 



T, . = u 



J'Dp 

167TW 



l/(2-/3 2 ) 



with 



D =Fpsm(ir/3 2 /2) 



r(/? 2 /4)r(i - /? 2 /2) 



-i 2 



r(i - /3 2 /4) 



(29) 



(30) 



Given the fact that we used quantitative estimates for 
the one-dimensional correlation functions and not just 
asymptotic estimates (^9|) should even be able to give 
semi-quantitatively the T c if J' can be determined by 
an independent method. Much more important however, 
is the field dependence of the T c . Indeed since the expo- 
nents depend on the field in a non trivial way one can ex- 
pect a non trivial magnetic field dependence of the three 
dimensional transition temperature. Close to h — an 
analytical solution can be obtained. We have 



d\nT c /J' 
dH 



3 df3 2 
~A~dH 



l.W + \n(v/T c a)} 



1 d\j\v 



4 dH 



(31) 



It is clear that at small enough J' this expression becomes 
positive and thus magnetic field increases the transition 
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temperature. A numerical estimate of T c is given in Fig 
ure 71. 
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FIG. 7. Three dimensional transition temperature T c 
as a function of the applied magnetic field h (i.e. for 
hd < h < h C 2). The full line is the complete solution whereas 
for the dotted curve the prefactor Fp has been fixed to its zero 
field value. Note the minimum at h — 

The camel-type structure, instead of the naively ex- 
pected dromedary one, comes from the competition be- 
tween the decrease of the exponent when one moves closer 
to h c i or h C 2 (which leads to an increase of the T c ), with 
the fact that the excitation velocity u decreases at the 
same time. 

We can also obtain the local susceptibility close to T c . 
In the mean field approximation it is given by 



dqi% 



(2tt) 2 J [x- 1 - 2 J[ - 2J' 2 ] + 4J[ sin 2 + AJ' 2 sin 2 q 2 /2 



When r = Xo^O, 0) - 2J[ - 2J' 2 < J' one can linearize the sinuses and calculate the integral 



In the limit of zero frequency we get 

lim xLH 

ui—tO LU 



■. tan 



o -i 



Mexo 1 ~ 2J( - 2J> 2 



dx{^{\ - /3 2 /4 + ix) - i/j(P 2 /A + ix)} 



1 - (T c /T) 



2-f3 2 



r(/3 2 /4~iz)r(l-/3 2 /4) 
r(l-/3 2 /4-ia;)r(/3 :4 /4) 



(32) 



(33) 



(34) 



Close to T c ( p4| ) reduces to a mean field divergence 
1/(T - Tc) 1 / 2 . This leads to a similar divergence in the 
relaxation rate 1/T±. 

B. Low density limit 

When the magnetic field is close to h c \ or h c2 the above 
mean field approach on the single chain Luttinger liquid 
cannot be used. Indeed the energy of the interchain cou- 



pling becomes larger than the intra chain energy scale, 
and interchain coupling should be treated from the start. 
Fortunately the problem is still solvable since the num- 
ber of excitations above the fully polarized ground state 
become very small. Let us focus on h ~ h c ± , the solution 
for h ~ h C 2 can be deduced by symmetry. 

Another useful way of viewing this problem, specially 
useful when we deal with the low density limit, is given 
by using the boson representation of spins, instead of 
the standard Jordan- Wigner fermionic one (the fermions 
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have to carry a string) . The spins can be represented by 
hard core bosons. The presence of a boson denotes a spin 
up state or in the original ladder a triplet on the rung. 
The hard core constraint ensures that one has only two 
states (full or empty) on each site to get a faithful repre- 
sentation of the spin 1/2. The problem thus reduces to 
a problem of hard core bosons with in chain interactions 
(due to the S Z S Z term). The interchain coupling Hamil- 
tonian is thus just the kinetic energy interchain hopping 
term of these bosons 



4 



E 

i,(a,0) 



b\, a b 



i,0 



(35) 



Since the bosons are very diluted it is essentially exact 
to neglect the interactions between them but for the hard 
core constraint, as indicated by the fact that the Lut- 
tinger liquid parameter for a single chain goes to K = 1 
close to h c i . One has thus to solve the problem of a three 
dimensional gas of hard core bosons with the simple ki- 
netic energy 



Jii 



E( 6 LA+^+ h - c -) 



(36) 



J± 



i,(a,/3) 

To go back to the standard negative hopping one must 
perform along the chain the gauge transformation 



(-1)** 



(37) 



Since the density is low (|36|) can be reduced to the con- 
tinuum limit of bosons with the kinetic energy 



k 2 



k 2 



(38) 



The ordering transition thus reduces to the well known 
problem of the iBpse condensation transition for the di- 
luted boson gasci The three dimensional ordered phase 
corresponds to the superfluid one, whereas the high tem- 
perature phase is the normal fluid of bosons. Most phys- 
ical quantities relevant for the ladder problem can im- 
mediately be borrowed from the vast knowledge existing 
for the diluted boson gas. We will not dwell on all quan- 
tities that can be computed but simply give here a few 
examples. 

The critical temperature is knowu_in the limit of low 
density and is given by the equatiorO 



A 

T 



h 
T 



(4tt) 3 / : 



r C(3/2)T 3/2 = 



(39) 



where to is a simple number (scattering matrix for an 
infinite hard core potential). This leads to a critical tem- 
perature varying as 



T c oc (h — h c i) 



2/3 



(40) 



and gives back immediately the mean field critical expo- 
nents. Of interest is also the total density p, i.e. the 
magnetization of the ladder. It is given by 



1 



P : 



(4tt) 3 / 2 



C(3/2)T 



3/2 



T < T c (41) 
T > T c (42) 



Thus two effects occurs. First the magnetization is non 
monotonous in temperature and increases by a factor of 
two between T — T c and T — 0. Second, at criticality 
h = h c i the magnetization grows as T 3 / 2 at very low 
temperatures (which can be readily seen by computing 
J d d qnB(e(q))- This is a different temperature depen- 
dence than for independent ladders. In that case it would 
be given by the excitations of the one-dimensional theory, 
which arc fcrmionic in nature and have a one dimensional 
density of state, leading to 



m cx / deN ld (e)n F (e) cx T 1/2 



(43) 



Looking at the temperature dependence (for very low 
temperatures) of the the magnetization at criticality 
should thus provide information on the interladdcr cou- 
pling. Of course at higher temperatures one always re- 
covers the independent ladders behavior ([43]). 

Another important quantity is the NMR relaxation 
rate. The (S + S~) correlation function is here simply 
given by the bosonic single particle Green's function 



(S+(r,r)S- (0,0)) = (-ir(c(r,r) C t(0,0)) 



(44) 



The (— l) r factor coming from the transformation (|37| ) 
implies that the low energy (massless) part of the spec- 
trum which for the spins is around q ~ tt (see Figure ||)) 
is given by the small oj small q Green's function for the 
bosons. The single particle Green'sJunction can also be 
computed in the limit of low densityCJ and is given in the 
condensed low temperature phase by 



G(u} n ,q) 



iuj n + k 2 + A 
uj 2 + k 4 + 2Ak 2 



(45) 



Using the standard formula for the NMR relaxation rate 
(O) and (E|) one obtains 



1 T 
-—- oc — == 

Ti VA 



(46) 



where A is given by (|39|) giving thus a relaxation rate pro- 
portional to T at low temperature and diverging close to 
T c as 1/(T C - T)V2. 

The above results could apply to the three dimen- 
sional phasejOf Cu 2 (C5H 12 N2)2Cl4, which is a strong cou- 
pling laddeili3 where this theory is directly applicable. 
Whether or not the transition experimentally observed 
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is due to the mechanism presented here is still an open 
question . and . other mechanisms of instability have been 
proposed 2 ^rT Various experiments can be performed to 
elucidate this point. First since the 3D transition de- 
scribed here is simply an ordering in the direction of the 
spins in the XY plane its impact on the global global 
magnetization is very weak, as seems to be the case ex- 
perimentally. Note however that it does change the tem- 
perature dependence of the magnetization at criticality 
and below T c . Other interesting experiments could be a 
fit of the h— T c relation (40) close to h c i and more gen- 
erally the camel-like shape of the phase diagram. Local 
probes like NMR or neutrons should be perfectly suited 
to study this transition. NMR could provide a way to 
map the phase boundary (by looking at the divergence 
of l/Tx ~ 1/|T C - Tl 1 / 2 ). The l/T x ~ T law in the 
low temperature phase could also provide conclusive evi- 
dence. Finally since we know that the transition is in the 
normal-superfluid transition universality class, one could 
also try to compare the thermodynamic singularities. 



V. CONCLUSION 

We have examined in this work the properties of lad- 
ders under magnetic field, and focussed on the gapless 
phase occuring between h c \ < h < h C 2- For a single 
ladder we computed quantitatively the correlation func- 
tions as a functions of the magnetic field. The correlation 
functions in the ladders have an identical structure both 
for the weak coupling ladder J± <C J\\ and the strong 
coupling one. As in weak coupling an interesting feature 
in the spin correlation function of the ladder compared 
to a single chain is the appearance for h ~ h c i of a low 
energy mode only close to q ~ for the (S Z S Z ) corre- 
lation function and close to q ~ n for the (S+S-) one. 
A single chain would have had both q ~ and q ~ tt 
mode at low energy. This prediction should be testable 
in neutrons experiments. Close to the lower critical field 
h c i where the, gap is destroyed we recovered the univer- 
sal exponentlia for the divergence of the NMR relaxation 
rate 1/Ti - 1/T 1 / 2 . 

These results served as a basis to analyze the nature 
of the phase transition in a system of three dimension- 
ally coupled ladders. This problem falls into the category 
of Bose condensation of hard core bosons, which allows 
to obtain many properties of the phase diagram and the 
ordered phase. The variation of T c with the field has a 
local minimum at h — (h c i + h C 2)/2 leading to an un- 
usual camel-like shape for the phase diagram. Close to 
h c i and h C 2 the transition is similar to the one of a di- 
luted Bose gas with T c ~ (h — ft, cl ) 2 / 3 . The temperature 
dependence of the S s magnetization goes at criticality 
h = h e i from a T 1 / 2 behavior for independent ladders 
to a T 3 / 2 (when interladder coupling intervenes) at low 
temperatures. The NMR rate diverges close to the tran- 
sition as 1/Ti ~ \T-T C \~ 1/2 and behaves as l/T x ~ T at 



low temperatures. These quite distinct features could be 
used to check whether this transition is the one occuring 
in the experimental system Cu2(C5Hi 2 N 2 )2Cl4. 



ACKNOWLEDGMENTS 

We thank L. P. Levy for many inspiring discussions. 
A. T. is grateful to University of Paris-Sud, where part 
of this work was completed, for support and hospital- 
ity. T.G. thanks the I.T.P. (Santa Barbara), where the 
final stage of this work was completed, for support and 
hospitality. This research was supported in part by the 
National Science Foundation under grant PHY94-07194. 



APPENDIX A: WEAK VS STRONG COUPLING 
LADDER 

We show in this appendix the connection between the 
correlation functions for the weak coupling ladder and 
the strong coupling one. The use of-the simple luttingcr 
liquid expressions due to Haldanecd allows for a more 



transparent derivation than the one given in Ref. 22 



For weak J± one introduces two boson fields (one for 
each leg) and it is more convenient to use the symmetric 
and antisymmetric combinations. 



± ■ 



n,2 



V2 



(Al) 



With the usual representations of the spin operators in 
terms of the (j) and 9. Since the magnetic field couples 
only to the symmetric field, it can only remove this gap 
and the antisymmetric field remains massive even above 
h c i- To write the correlation functions for the spins in 
terms of the bosonic fields one use the standard represen- 
tation of spin in terms of fermion operators and then the 
expression of this fermions in terms of the bosonic ones 
for a Luttinger liquid. Let us start with the S z operator, 
connected with the density of associated fermions. 

We use the decomposition of the density or fermion 
operator in a Luttigner Liquid. 



P = Po- -V a 

7T 



in(Qx+4 



(A2) 



which contains all harmonics 2kp, Akp etc. of the fermion 
density. Traditionally one only retains the lowest (most 
singular) harmonic, which leads to the standard expres- 
sion (|^). However here, since here the 2k f component is 
massive due to the presence of the antisymmetric mode, 
it makes it necessary to retain the next harmonic. This 
leads to the spin operators (e.g. for spins on chain 1) 

S z = V0 S + e l( 3 a; + v / 2(0 s +0a)) _|_ e i2Qx+2V2(4>s+<j>a) _| 

(A3) 



9 



Since the field <p a remains massive even above h c i , all cor- 
relation functions containing it decay exponentially, and 
can be neglected at large distance (or time). Thus no 
Q component appears in the correlation function for the 
laddeiO, in marked contrast to the frustrated or dimer- 
ized single chain. On the other hand the 2Q term con- 
tains 2y/2(f> a . Although this term is superficially massive, 
it can be combined with a cos(2\/2</>a) term existing in 
the Hamiltonian for the ladder, giving rise to the opera- 
tor 



Ce 



i(2Qx+2V2{<p e )) 



(A4) 



where c is a mere constant. This operator containing only 
the symmetric field is massless. The long wavelength de- 
cay of the correlation function in the weak coupling lad- 
der is thus given by 



(S z (r)S z (0)) 



1 



C 2 cos(2Q) 



4 K 



(A5) 



For the weak coupling ladder since K > 1/2 the cos(2Q) 
term is subdominant and can be safely dropped. To 
make the connection with the case of strong .coupling, 
where one can have K < 1/2, it must be kept£3. Since 
Q = 7r(l — 2m) it is easy to see that ( A.5) has exactly the 
same form than the expression derived for the strongly 
coupled ladder (|lj), showing that the two limits are 
smoothly connected. 

Similar results can be obtained for the higher harmon- 
ics 2nQ. For the transverse magnetization correlation 
one get in a similar way for the weak coupling ladder 

S+ = e te [{-I) 1 + cos(20) + (-I) 1 cos(40)] (A6) 

where the cos(4</>) term comes again from the higher 
harmonics. As for the S z component, the cos(2</>) re- 
mains massive due to the gap in the <pa field, whereas 
the cos(40) can again be combined with terms in the 
Hamiltonian to give a massless term. The final result is 



(s»s-(o)> = {-iy\- 



1/iK 



cos(7r(l + 2m)r) ( - 



1/4K+4K 



(7) 



Thus the expression ( J14| ) for the strongly coupled lad- 
der is again similar to this one. 

Thus weak and strong coupling ladders are smoothly 
connected. The crucial reason is that the gap in the anti- 
symmetric degrees of freedom which exists already in the 
weak coupling ladder is equivalent to the neglect of the 
two excited states of the triplet performed for the strong 
coupling. In a system without such an antisymmetric 
gap (such as a dimerized chain) this smooth continuity 
would not hold and the weak and strong coupling corre- 
lation functions would be radically different. 
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